7V7  9-^f 


MISCELLANEOUS  PAPER  NO.  1-808  v - 

WATER  SHOCK  WAVES  FROM 
ABOVE-WATER  EXPLOSIONS 

by 

J.  M.  Pinkston,  Jr. 

Akim  Sakurai 


March  1966 


i 


R*produtcrt  by 

NATIONAL  TECHNICAL  ' 
INFORMATION  SERVICE  ' 

U  $  Ofpuftmenf  of  CommtfCt  I 

SiMingfu-lil  VA  72151 


U.  S.  Army  Engineer  Waterways  Experiment  Station 

CORPS  OF  ENGINEERS 


Vicksburg,  Mississippi 


t 


O 

t *:J 


pa 

•**] 

N<J 

O 

he.1 


j 

) 

c*-j 

Cl 

b’—i 

1 

L-  j 

t-  * 

*«• . , 

t.1  - 

M.] 

IT  i 
L:  .> 

« 

C.T 

.....  J 
.! 

f'  ;•  ;1 

»;  ;n 

f-  -J 

M 

t—  4 
.  - 

Utr  -j 


UBRARY 

US  ARMY  ENGINEER  WATERWAYS  EXPER  ME: IT  STATION 

yiCKSouiu;  iaississjs’pj 


<* 


MISCELLANEOUS  PAPER  NO.  1-808 

WATER  SHOCK  WAVES  FROM 
ABOVE-WATER  EXPLOSIONS 

by 

J.  M.  Pinkston,  Jr. 

Akira  Sakurai 


March  1966 


U.  S.  Army  Engineer  Waterways  Experiment  Station 
CORPS  OF  ENGINEERS 


••mv.mmo  viOKaauaa,  mi«b. 


Vicksburg,  Mississippi 


SanrfSiTVr  Jtg  ;v  S-  'F^T^i 


mMr 


r/f  7 


-.  5  J. 


j  *  Jf  Av-*  v 


|  '  Foreword 

5 

|  xhis  paper  was  prepared  for  presentation  at  the  1S)66  Army  Science 
9.  fjcrcnce  at  West  Point,  New  York.  It  was  approved  for  presentation  and 

I  Plication  by  the  Office,  Chief  of  Engineers.  The  paper  was  prepared  by 
J.  M.  Pinkston,  Jr.,  and  Dr.  Akira  Sakurai,  under  the  general  super¬ 
men  of  Mr.  G.  L.  Arbuthnot,  Jr.,  Acting  Chief,  Nuclear  Weapons  Effects 
>:*lslon,  and  Mr.  J.  N.  Strange,  Chief,  Engineering  Research  Branch, 

'  S.  Army  Engineer  Waterways  Experiment  Station. 

Director  of  the  Waterways  Experiment  Station  during  the  preparation 
,-f  this  paper  was  Col.  John  R.  Oswalt,  Jr.,  CE.  Technical  Director  was 
Jj.  J.  B.  Tiffany! 


I 

[?!i 

d 


tr  r,  n  n  a 


*  \£%i r^flrrie^'  $3^*3'^: 


tt*9*SH?Z)£* 


/vf 

'1 


PINKSTON. AM)  SAKURAI 

TITLE:  Water  Shock  Waves  from  Above-Water  Explosions  ! 

J.  M.  PINKSTON,  JR.,  AND  AKIRA  SAKURAI 
U.  S.  Army  Engineer  Waterways  Experiment  Station 
;  Vicksburg,  Mississippi 

I 

ABSTRACT:  j 

i 

In  an  effort  to  determine  the  amount  of  energy  trans¬ 
mitted  into  the  water  when  an  explosion  occurs  above  a  water  sur¬ 
face,  both  experimental  and  theoretical  investigations  were  con¬ 
ducted.  The  experimental  study  involved  detonation  of  21-lb 
spheres  of  TNT  above  a  water  surface  and  measuring  induced  water 
shock  at  spatial  positions  underwater.  To  fully  incorporate  the 
nonlinear  characteristics  of  the  airblast  into  the  study,  free- 
air  pressure-time  data  were  used  as  the  water  surface  loading 
mechanism  and  the  investigation  was  divided  into  two  phases— in 
i Phase  I  the  generating  source  was  at  such  a  height  that  the  in- 
■duced  water  shock  could  be  adequately  described  by  acoustic 
theory,  while  in  Phase  II  the  generating  source  was  nearer  the 
| surface  and  the  resulting  disturbance  field  was  of  a  finite  am- 
1 plitude  that  had  to  be  considered  to  accurately  determine 
junderwater  pressure. 

j  Peak  pressures  and  pressure-time  histories  derived  dur¬ 

ing  P.iase  I  (acoustic  theory)  compare  very  well  with  experimental 
'.results  for  those  cases  wherein  the  theory  is  valid.  When  the 
generating  source  nears  the  surface  however,  the  theory  becomes 
inadequate  (primarily  because  of  a  rigid  interface  assumption) 
|and  theoretical  results  are  adjusted  by  a  correction  factor  (0) 
that  modifies  the  pressure  values  as  0p  .  Appropriate  0  val¬ 
ues  for  each  charge  position  were  obtained  by  comparing  theoreti- 
i cal  and  experimental  results  for  each  charge  position. 

!  Phase  II  results  (the  effect  of  the  finite  disturbance) 

. indicated  that  the  disturbance  is  not  as  important  as  originally 
I  estimated,  and  except  near  the  zero  point,  the  acoustic  theory 
can  offer  a  complete  description  of  the  pressure  phenomena  in 
water.  However,  near  the  zero  point,  the  finite  disturbance  ef¬ 
fect  is  essential  to  an  accurate  determination  of  the  underwatc.: 
(pressure. 
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PINKSTON  AND  SAKURAI 

'developed  during  the  investigation  would  be  ' in  such  a  form  that  the"  “ 
numerical  procedure  to  determine  the  pressure-time  history  at  spatial 
positions  underwater  could  be  easily  performed. 

■  The  theoretical  investigation  was  divided  into  two  phases. 

•In  Riase  I  it  was  assumed  that  the  generating  source  was  at  such  a 
•height  above  the  surface  that  the  induced  water  shock  could  be  ade¬ 
quately  described  by  acoustic  theory,  while  Phase  II  covered  those 
'cases  in  which  the  generating  source  was  nearer  the  surface  and  the 
i  resulting  disturbance  field  in  water  was  of  a  finite  amplitude  that 
■'~;had  to  be  considered  in  determining  underwater  pressure. 

|  PHASE  I  -  ACOUSTIC  THEORY 

4 

».  {  #  » •  ••  -- 

.  ;  In  this  phase  of  the  investigation  the  problem  can  be  re- 

\  duced  to  that  of  finding  the  pressure  field  (p)  in  the  water  from 
•the  wave  equation,  subject  to  the  boundary  condition  given  by  the 
•pressure  loading  pQ  on  the  surface.  Assuming  Pq  as 


P0(r>t)  = 


'a0(t)-  +  r^^t)  +  r*a2(t)  +...,  for  R  >  r  >  0 


(1) 


i 


for  r  >  R 


1 


(see  Figure  l) 


>  where  &o(t),  ...,  are  functions  of  time  (t),  the  general  solution 
|  for  p  can  be  simplified  into  an  expression  given  by  a  series  of 
6in6le  integrals  including  aQ(t),  a^t),  ...,  in  each  integral  as 

'  ;  *  p0(r,t  -  i}'  - 


ri  r 2 


A - 


_t2(r,z,t) 

L  /  4^r»z,t»t^an(t^dt‘'  (2) 

n=0  t1(r,z,t) 


;  ’-'here  H  is  the  Heaviside  function  and  /  are  known  functions  "of" 
|  r  *  z  ,  t  ,  and  t’  (z  =  gage  depth).'  Thenlimits  t.  and  t2  are 

•  functions  that  are  determined  for  each  gage  location  (r.z)  and'  * 
i  time  t  . 

!  ao(t)>  •••>  functions  were  determined  using  both 

•  experimental  free-air  pressure-time  data  and  the  properties  of  a 

•  reflected  airblast  wave  at  a  rigid  interface.  Although  uractically 
:  Any  number  of  these  functions  can  be  determined,  indications  are1, 

!  Ji'At  when  compared  with  experimental  data  sufficient  accuracy  is 
utained  using  only  the  second  approximation,  i.e.  retaining  just 
‘  and  al(t)  terras.  These  simplifications  make  it  easier 

.  determine  the  pressure-time  history  at  a  given  positic; 
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'underwater  by  numerical  integration  over  the  desired  time  interval.""] 


peak  pressure  values  were  determined  directly  from  Equation  2  using 
a  limiting  process.  All  equations  were  programmed  for  computer 
processing,  and  the  results  obtained  are  compared  with  experimental 
data  in  Figures  2  through  7. 

Since  the  acoustic  theory  becomes  inadequate  when  the 
generating  source  nears  the  water  surface  (primarily  because  of  the 
rigid  interface  assumption  used  to  determine  the  aQ(t) ,  ,  func¬ 

tions),  theoretical  values  given  by  the  acoustic  theory  are  not  in 
good  agreement  with  experimental  results;  however,  these  theoretical 
values  can  be  partially  adjusted  using  a  correction  factor  (0)  that 
modifies  the  pressure  values  as  0p  .  Appropriate  0  values  (which 
should  be  in  the  range  of  0  <  0  <  l)  for  each  charge  position  were 
obtained  through  a  formal  comparison  of  theoretical  and  experimental 
results  as  shown  in  the  following  paragraphs. 

Peak  pressure  values  for  the  case  when  the  charge  is 
directly  above  the  gages  (Xj.  =  0)  are  plotted  (solid  lines)  in 
Figure  2.  The  upper  half  of  the  figure  illustrates  test  geometry 
where  Xjj  is  reduced  charge  height,  Xr  is  reduced  horizontal  dis~- 
tance  from  charge  surface  zero  to  gage  surface  zero,  and  z  is  the 
gage  depth  in  feet.  The  two  curves  for  each  X«  correspond  to  0 
values  of  0.8  and  1.0.  For  reduced  charge  heights  of  4x  and  5X  , 
theoretical  and  experimental  results  are  in  excellent  agreement, 
but  as  suspected,  the  theory  becomes  inadequate  as  the  charge  nears 
the  water  surface.  Figure  2  establishes  the  range  and  validity  of 
the  present’  theory  and  provides  an  estimate  of  the  appropriate  0 
value  for  each  charge  position.  (A  0  value  of  0.6  when  Xjj  =  1 


and  0.7  when  XH  =  2  or  3.) 


Figure  3  is  a  plot  of  arrival  time  as  a  function  of  gage 
depth  and  indicates  that  the  shock  wave  in  water  is  well  documented 
by  acoustic  theory.  Peak  pressure  calculations  for  spatial  posi¬ 
tions  underwater  were  equally  successful  when  compared  with  exper¬ 
imental  data  as  shown  in  Figures  4  and  5  for  the  illustrated  test 
conditions . 

Experimental  pressure-time  results  are  compared  with  theory 
in  Figure  6  for  the  case  when  the  charge  is  directly  above  the  gages  I 
(Xr  =  0).  The  results  compare  very  well,  particularly  for  the 
deeper  gages.  Pressure-time  histories  at  spatial  positions  under¬ 
water  are  contained  in  Figure  7.  Except  for  an  apparent  difference 
in  shock  wave  decay  rates  for  the  deeper  gages,  theoretical  and 
experimental  results  are  in  good  agreement  and  further  improvements 


can  be  expected  by  retaining  additional  terms,  that  is,  a. 


in  the  basic  equation. 


PHASE  II  -  THE  FINITE  DISTURBANCE  EFFECT  IN  WATER 


As  stated  earlier,  when  the  generating  source  nears  the 
water  surface  the  acoustic  theory  is  no  longer  valid  and  the  effect 
of  the  finite  disturbance  must  be  considered.  Though  the  distur¬ 
bance  is  of  a  finite  amplitude,  it  is  reasonable  to  assume  that  the 
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'disturbance' field  is  restricted  to  a  small  region  on  the  water  surface, 
figure  8  illustrates  this  particular  case  which  is  characterized  by  j 
:  interaction  of  the  airblast  and  water  shock  fronts  at  point  c  on 

’the  water  surface.  When  this  occurs,  the  water  is  compressed,  its 
1  front  is  a  genuine  shock  wave ,  and  a  nonlinear  analytical  approach 
i  required  to  determine  the  flow  field.  Fortunately  the  flow 
■  field  near  the  intersection  (point  c)  acts  independently,  and  its 
1  exact  solution  can  be  determined  through  the  Rankine-Hugoniot  rela- 
i  tions  of  both  the  air  and  water  shock  fronts  and  knowing  the  condi¬ 
tions  that  exist  at  the  interface.  These  relations  are  reduced  to 
|  a  system  of  transcendental  equations  to  determine  the  angles  of 
'  intersection,  in  addition  to  pressures,  velocities,  and  densities 
;  in  both  the  air  and  water.  The  equations  were  solved  numerically, 

;  and  the  water  shock  pressure  values  obtained  for  various  reduced  _ ■ 

I  charge  heights  are  shown  in  Figure  9*  (These  values  correspond  to 
i  Pq  values  when  r  =  R  in  Equation  1.) 

1  Table  1  compares  p0  values  for  r  =  R  =  0  from  Phase  I 

|  with  corresponding  values  from  Phase  II.  Since  Phase  I  results  are 

i  based  on  "the  rigid  interface  assumption,  the  difference  in  the  two . 

[-values-  indicates  what  effect  the  compressibility  of  water  has  on 
its  reflection  characteristics. 


TABLE  1 


X  PRESSURES  AT  r  =  R  »  0  FOR  VARIOUS 


VALUES 


0.1 

0.5 

1.0 

2.0 

3.0 

”  4.o 

67,390 

16,212 

6,566 

1,872 

594 

291 

63,763 

15,697 

6,432 

1,852 

591 

291 

Phase  I 
Riase  II- 


Results  show  a  decrease  in  pressure  that  is  more  pronounced 
!  near  the  water  surface.  Though  the  difference  is  rather  small  and 
j  perhaps  insignificant,  except  for  the  Xjj  =  0.1  case,  these  data 
.  do  suggest  that  the  effect  of  the  finite  disturbance  in  water  is 
1  not  as  important  as  originally  estimated. 

•  Using  the  solution  near  the  point  c  ,  a  successive  ap- 
!  proximation  method  to  determine  the  entire  flow  field  in  air  and 

|  water  was  developed  and  its  first  approximation  studied  in  detail. 

,  Figure  10  shows  a  plot  of  pressure  values  from  this  first  approx- 
|  imation  for  the  Xjj  =  0.5  »  X~  =  0  case  as  compared  with  available 
j  experimental  data.  Again,  indications  are  that  the  acoustic  theory 
;  can  offer  a  complete  description  of  the  pressure  phenomena  in  the 
j  water,  with  the  exception  of  the  near  zero  point,  if  the  appropriate 
j  0  value  is  selected.  In  the  vicinity  of  the  zero  point,  however, 

:  the  finite  disturbance  effect  in  water  is  essential  to  an  accurate 

•  determination  of  the  underwater  pressure. 
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Figure  5,  'Peak  pressure  at  various  gage  locations  for  a 
reduced  horizontal  distance  (X_)  of  4 
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